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Abstract 

The existence and nonexistence of global in time solutions is studied for 
a class of equations generalizing the chemotaxis model of Keller and 
Segel. These equations involve Levy diffusion operators and general 
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1 Introduction 

We consider in this paper the following nonlinear nonlocal evolution equation 
generalizing the well known Keller-Segel model of chemotaxis 

dtu + {-A)''/^u + V ■ {uB{u)) = 0, (1.1) 

for (x, t) € M"^ xM"'', where the anomalous diffusion is modeled by a fractional 
power of the Laplacian, a € (1,2), and the linear (vector) operator B is 
defined (formally) as 

5(n) = V((-A)-^/\). (1.2) 

For /3 € (1, d], and d >2, one may express the nonlocal nonlinearity in p.ip 
using convolution operators since 

B{u)ix) = sd,p / -r—TdZ^u{y) dy, (1.3) 
jRd |x — y\ 



1 



with some Sd,/3 > 0, and the assumption /3 > 1 is needed for the convergence 
of this integral. 

Of course, the choice a = 2, /? = 2 in (ll.ip - (ll.3p corresponds to the usual 
Keller-Segel system studied mainly in space dimensions d = 1, 2, 3. It is 
well known that if /? = 2, the one-dimensional system (ll.ip - ()1.3p possesses 
global in time solutions not only in the case of classical Brownian diffusion 
a = 2 but also in the fractional diffusion case 1 < a < 2 as was shown 
by C. Escudero in [17J. On the other hand, there are many results on the 
nonexistence of global in time solutions with "large" initial data if d > 2 and 
a = 2, /3 = 2, see, e.g., [H [H H [H] . Even if d > 2, a = 2 and 1 < /3 < d, 
there are results on the blow up of solutions with suitably chosen initial 
data, see |10] (caution: the notation in [10, Prop. 4.2] differs from that in 
the present paper). Let us finally recall that, in the limit case a = 2, (3 = d, 
mass M = J^aUo{x)dx is the critical parameter for the blow up, see 
Prop. 4.1] and [H]. 

The usual method of proving the nonexistence of global in time nonneg- 
ative and nontrivial solutions, used in the abovementioned papers, consists 
in the study the evolution of the second moment of a solution 'W2{t) = 
J^^\x\'^u{x,t) dx and to show (via suitable differential inequalities) that 
W2{t) vanishes for some t > 0. The second moment of a typical solution 
to an evolution equation with fractional Laplacian cannot be finite, see e.g. 
[13j . Hence, our goal in this paper is to generalize the classical virial method 
and to show the blow up of solution system (ll.ip - ()1.3p by studying moments 
of lower order 7 G (1,2) 



After this paper was completed, we discovered a recent preprint [23] 
where the authors show the blow up of solution to system (jl.ip - (jl.3p with 
fractional diffusion in the particular case d = 2 and j3 = 2. Our argument 
is different than that in [23j, shorter, seems to be more direct, and applies 
in more general situations. Moreover, we are able to formulate a simple 
condition on the initial data which leads to the blow up in a finite time of 




(1.4) 
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the corresponding solution. 

Notation. The L^-norm of a Lebesgue measurable, real- valued function 
V defined on R'^ is denoted by ||f The constants (always independent of 
X, t) will be denoted by the same letter C, even if they may vary from line 
to line. Sometimes, we write, e.g., C = C(*) when we want to emphasize 
the dependence of C on a parameter "*" . 

2 Main results 

The crucial role in the approach in this paper is played by the following 
scaling property of system (jl.ip - (jl.3p 

M^(x, t) = A°+^~2u(Ax, A"t) for ah A > 0, (2.1) 

in the sense that if n is a solution to ()l.ip - (ll.3p . then is so. In par- 
ticular, in our construction of solutions to ()l.ip - (ll.3p we use the fact that 
the usual norm of the Lebesgue space L"^/^"^^^^-* (M'^) is invariant under the 
transformation uq{x) i— > \'^^I^~'^uq{\x) for every A > 0. 

Theorem 2.1. Assume that d>2, a ^ (1)2], and /5 € {l,d\. Let 

( d 2d ] 
max < — -, — - > < p < d. 

\a + P-2' d + (3-l\ 



i) For every uq € LP(R ) there exists T = T(||no||p) and the unique local 
in time mild solution u G C([0, T], LP(IR'^)) of system il.l \) -U.3 \) with 
uq as the initial condition. 

ii) There is e > such that for every uq G L'^/^"+^~^^(M'^) satisfying 

\\uo\\d/{a+l3-2) < ^ (2-2) 

there exists a global in time mild solution u G C([0, oo), L^'(M'^)) of 
system !il.l\) - FT^) with uq as the initial condition. 

Moreover, if uq{x) > 0, then the solution u in either i) or ii) above is non- 
negative. Finally, if G L^(M'^), then the corresponding solution conserves 
mass 

/ u{x,t)dx= / UQ{x)dx = M. (2.3) 
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Recall that a > 1 is a usual assumption ([10, Th. 2.2], [HI [9]) which 
permits us to control locally the nonlinearity in (ll.ip - (ll.3p by the linear 
term. 

The results stated above can be easily generalized for equations with 
general Levy diffusion operators considered in [8, i9j but we do not pursue 
this question here. We refer the reader to the above mentioned papers, as 
well as [10], for physical motivations to study such equations. On the other 
hand, motivations stemming from probability theory {propagation of chaos 
property for interacting particle systems) can be found in, e.g., [7\. 

Remark 2.2. For a € (1,2) and (3 > I satisfying a + /3 > d + 2, it can 
be shown that local in time solutions can be continued to the global in 
time ones, see \10\ Th. 3.2]. However, in [10], another approach (via weak 
solutions) is used to construct solutions of system p.ip - (|1.3|) . 

The proof of Theorem 12.11 on local and global solutions to p.ip - (|1.3p 
follows a more or less standard reasoning which we sketch in Section [3l Our 
main goal, however, is to prove the finite time blow up of solutions to the 
nonlocal system (ll.ip ~ ()1.3p . 

Theorem 2.3. Assume that d> 2. The solution of ( fj.i| )-i f73j] with a non- 
negative and nonzero initial condition u{x,0) = Uo{x) blows up in a finite 
time in each of the following cases: 

i) (large mass) for a = 2, (3 = d, uq & L^{M.'^, (1 + |xp) dx), and if 



the threshold value of M is Sir if d = 2; 

a) (high concentration) for a € (1,2] and (3 € (l,(i] satisfying a + (3 < 
d + 2, uq & L^{R'^, (1 + l^l''') dx) for some 7 e (1, a), and if 



M = uo{x) dx > 2d/sd,i3, 
with the constant Sd^/3 defined in lil.3\) : in particular: §2,2 



~ 2n 



TT- SO that 




(2.4) 



for certain (sufficiently small) constant c > independent of uq. 
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The result stated in i) is essentially contained in [14J. The condition for 
blow up in the form ()2.4|) appeared already in [2j and [25], of course, for 
a = 2 only. Note that i) is a limit case of ii). Indeed, (12. 4|) written as 



becomes a condition on (sufficiently large) mass: 1 < cM, when (a + /?) y 



We have to emphasize that Theorem I2.3l ii contains a result which is 
new even for the classical parabolic-elliptic Keller-Segel model {i.e. equa- 
tions ()l.ip - (ll.3p with a = (3 = 2). Indeed, for d > 3, the conditions from 
part ii) of Theorem 12.31 guarantee the blow up in a finite time if the moment 
of order 7 of the initial condition is finite for some 7 S (1)2). All other 
known proofs of the blow up required just 7 = 2. An analogous result for 
d = 2 is stated below. 

Corollary 2.4. Assume that a = f3 = d = 2in equations lll.l\) - n~3\) . 

Suppose that there exists 7 € (1,2) such that uq G L^{M.'^,{1 + |x|'^)dx). 
There exists > such that if 



then the solution of lll.l\) - n~S\) with the initial condition u{x,0) = Uo{x) 
blows up in a finite time. 

It is conjectured that M2 = = Svr for every 7 € (1,2) and we 
expect that it can be shown by a careful analysis of constants appearing in 
inequalities ([OD - ijiTO]) and in (|il8]) - (fO0D in the proof of Theorem [231 
Here, one should recall that radially symmetric solutions to equations (jl.ip - 
()1.3p with a = P = d = 2 blow up in finite time under the assumption that 
their mass is larger than Svr and no moment condition imposed on the initial 
data is necessary, see [U Prop. 2.2, Th. 3.1] for a detailed presentation. 

Remark 2.5. Due to the translation invariance of problem (ll.ip - (ll.3p . the 
conditions on moments can be imposed on the quantity 




d+2-a-l3 



{d + 2). 




7' 
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where x = [j^d xuq{x) dx)/{J^d uq{x) Ax) is the center of mass of uq. 

Remark 2.6. Let us observe that the assumptions (j2.2|) and (|2.4p are in 
a sense complementary due to the following elementary inequality involving 
the LP-norms, mass M = J^^ u{x) dx, and the moment w.^ = \x\^u{x) dx 
of a nonnegative function u 

HI. >»/(-) . (2.5) 

To prove ()2.5p observe that > W Ld\ n (m u{x) dx, so that 



7 - JlRd\Bfl(0) 

/ u{x)dx = M- u{x)dx> M - R~^w^. 

Multiplying both sides of this inequality by R"^^^ ^) we get with C = lu^ ^ 

1 / \ i-i 



C\\u\\p = / uP{x)dx ] I dx] Vp 

KJr'^ J \Jbr{o) J 

> R'^i^A [ u{x)dx> R'^^^A M - R-^fw^ 



Taking the optimal R, i.e. R^ = Cwj/M, we get (12. 5p . 

Now, it is clear that if condition (12. 4p for blow up is satisfied for some 
7 € (1,2] and a suitable constant c, then ioT p = d/{a + (3 — 2) appearing in 
Theorem [2Tlii. we have ||n||p > CMM d+2-a-f3j\ pj ^ (j ^ condition 
(I2.2p for global existence is violated for sufficiently small e > 0. Vice versa, 
if (12.21) is satisfied, (12.41) cannot be true with small constants c. 



Remark 2.7. Note that one can prove in a similar way the inequality 

\\u\\Mv>CM\^-—j , (2.6) 

where the Morrey space norm ||u||Afp (< C'||n||p) is defined as 

sup ^) / \u{x)\dx, 

-R>0,a;oGR'* Jbr{xo) 

see [21 (15)] in the particular case 7 = 2. The scale of Morrey spaces 
is relevant to study mean field type problems related to the Keller-Segel 
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model, since the Morrey norms are particularly well suited to measure the 
(local) concentration of densities, see e.g. p]. 

We prove Theorem 12.31 and Corollary 12.41 by showing the extinction 
in a finite time of the function w{t) = J^a ^'rix)u{x,t) dx where ^p^yi^x) is 
smooth and behaves like \x\'^ with some 7 € (1,0?) for large see ()4.ip 
below. Note that if a < 2, we cannot expect the existence of higher or- 
der moments defined in (|1.4p with 7 > a. Indeed, even for the linear 
equation dtv + (— A)°/^t; = 0, the fundamental solution Pa{x, t) behaves like 
Pa{x,t) ~ (t"'/" + ^ , and therefore the moment (|1.4|) with 7 > a 

cannot be finite, see and references given there. Thus, we cannot apply 
the usual reasoning which involves an analysis of the evolution of the second 
moment 1112 of the solution because the integral defining 'W2 may diverge. 

We recall in Proposition 13.41 a result showing that the moment p.4|) is 
finite for a large class of initial conditions in the case of 7 < a. 

3 Existence of solutions 

We are going to construct solutions to system (ll.ip - p.3p via the following 
integral formulation 

u{t) = Sa{t)uo- [ V ■ Sa{t-T){u{T)Bu{T))dT, (3.1) 

Jo 

i.e., we consider mild solutions. Here Sait)uo = Pa{t) * "^^o is the solution to 
the linear Cauchy problem 

dtv + {-^)'''^v = {), v{x,{))=uo, (3.2) 

and Pa{x,t) is the fundamental solution of ()3.2p which can be represented 
via the Fourier transform Pa{^,t) = e"*'^'". In particular, 

Pa{x,t) =t-'^''^Pa{xr^/'^), 

where Pa is the inverse Fourier transform of e^'^'", see |18i Ch. 3] and |13] 
for more details. It is well known that for every a G (0, 2) the function Pa 
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is smooth, nonnegative, and satisfies the (optimal) estimates 

< Pa{x) < C{1 + |x|)-(°+'^) and \VPa{x)\ < C{1 + |2;|)-("+^+i) (3.3) 

for a constant C and all x G M'^. Hence, it follows immediately from the 
Young inequality for the convolution and from the self-similar form of the 
kernel Pa{x, t) that for every 1 < q < p < oo there exists C = C{p,q,a) > 
such that 

\\Sa{t)uo\\p < Cr-^^'^hWoWq (3-4) 

and 

||V5Q(t)uo||p < Ct~^(^~^)~^\\uo\\q (3.5) 
for every uq € L^(R'^) and ah t > 0. 

The construction of solution to (13. Ih in spaces is based on the follow- 
ing abstract result, see e.g. [22], [21]. 

Lemma 3.1. Let {X, \\ ■ be a Banach space and H : X x X ^ X 
a hounded bilinear form satisfying \\H{xi,X2)\\x ^ ??||2;i||a'I12;2||a' for all 
xi,X2 G X and a constant r] > 0. Then, if < e < l/(4r/) and if y X is 
such that \\y\\ < e, the equation u = y + H{u,u) has a solution in X such 
that \\u\\x < 2e. This solution is the only one in the ball B{0,2e). 

We skip an easy proof of this lemma which is a direct consequence of the 
Banach fixed point theorem. 

Sketch of the proof of Theorem \2.1[ The proof consists in constructing solu- 
tions to the "quadratic" equation (j3.ip using Lemma [STT] with y = Sa{t)uo 
and with the bilinear form 

H{u,v) = -[ V ■ Sa{t-T){u{T)Bv{T))dT. (3.6) 
Jo 

Local existence of solutions. It suffices to obtain estimates required by 
Lemma |3. II in the Banach space Xj, = C{[0,T], L^{M.'^)) supplemented with 
the usual norm = sup^gjo^p] ||u(t)||p. 

By inequality (|3.4p . we immediately obtain ||»S'q,(-)'Uo||;^'P < ||wo||p- 
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Combining the definition (II. 3p of the operator B with the Hardy-Little- 
wood-Sobolev inequahty we obtain 



B{u)\\g < Sd,p\\ I • 1"''+^-^ * u\l < C\\u\\j, (3.7) 
1 _ §^ = I 

P d. q- 



for every 1 < p < q < oo satisfying ^ — ^-p- = i. Moreover, inequahty p.Sp 



and the Holder inequahty lead to 

||V5„(t - T){uB{vmp < C{t - r)-^-^(^-F) \\uB{v)\\r 

<C{t-T) « °V'- p)\\u\\p\\B{v)\\g. 

where ^ = | + |- Hence, by inequalities (j3.7|) and (j3.8|) . there exists a con- 
stant C such that for all n, u E X^, the bilinear form (j3.6p satisfies 

< sup / (t-r) " "V- W||n(r)||p||7;(T)||pdr 

tG[0,T] Jo (3.9) 

1 d 1 1 ^ 



II II iO,^ II II ^T.^ 

In estimates (13. 9p , we have used the relations 7 = ^ + | = | — > and we 
assume that 

• p > d/{a + P — 2) in order to have 1/a + {d/a){l/r — 1/p) < 1; 

• p > 2d/ {d + P — 1) to guarantee that r > 1; 

• P < — 1) to be sure that r < p. 

Choosing a sufficiently small T > in (j3.9p we complete the proof of 
Theorem 12. H i by an application of Lemma l3.ll 

Global in time solutions. Here, the reasoning is completely analogous: 
using inequalities p.4p . (j3.5p . ()3.7p - ()3.9p we estimate the bilinear form (j3.6p 
in the Banach space 

=C7([0,oo),L'^/(°+^-2)(M'^)) 

d / 1 a + l3-2 \ 

n{u G C7((0,cx)),LP(M'^)) : supt^Vp ^ / < oo} 

t>o 

supplemented with the norm 

d f 1 Qf+;3-2 A 

= sup||'u(t)||rf/(„+^_2) +supt°Vp d 

t>0 t>0 
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We skip further details of this standard reasoning. 

Nonnegativity property. In order to prove that uq > imphes u{t) > 0, it 
suffices to study the function u~{x,t) = max{— i), 0} and to fohow the 
arguments either from [9, Prop. 3.1] or from Prop. 2] in order to show 
that u~{x,t) = 0. Here, we do not give a detailed presentation because the 
proof from p3i Lemma 2.7] can be rewritten in this more general case. 

Conservation of the integral. If uq € L^(^'^), one should repeat the fixed 
point argument from i) (or from ii)) in the space Xj. fl C([0, T], L^(M^)) (in 
AfPnC([0,oo),Li(M°')), resp.) in order to have ^(t) G ^^(R"') for alH G [0,r] 
{t G (0, oo), resp.). Next, it suffices to integrate over the both sides of 
equation (|3.1|) . Using the following consequences of the Fubini theorem 



/ 



Pa{t) * uq{x) dx = I UQ{x)dx 
JW'- Jr'' 

and 

/ Vpa{t) *v{x)dx = for every v G L^{R'^), 
we conclude that J^d u{x, t) dx = J^d uo{x) dx for every i G [0, T]. □ 

Remark 3.2. The reasoning from the proof of Theorem 12.11 follows the lines 
of the usual proof of the local in time existence (as well as the global in time 
existence for small initial conditions) of solutions to system (jl.ip - fjl.Sp with 
a = 2. Since that argument is based on an integral representation analogous 
to that in (j3.ip and on counterparts of decay estimates from (j3.4p - (j3.5p . it 
can be easily adopted to the more general case of a G (1,2]. Examples 
of such a reasoning applied to various semilinear models and realized in 
miscellaneous Banach spaces can be found in [3l [5l [6l [TTl [201 IZH Ell [23] . 

Next, we recall weighted estimates of solutions to the linear Cauchy 
problem (13.20 which have been proved in, e.g., [13]. In the following, we use 
the weighted spaces 

L^{R'^) = {ve L°°(]R'^) : H^^Hl^o = esssup^gKd(l + \x\f\vix)\ < oo} 

for a fixed > 0. 
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Lemma 3.3. ([13, Lemma 3.1]) Assume that vq G -L^^(M'^). There exists 
C > independent of vq and t such that 

\\Sait)vo\\L^^^ < C(l +t)||vo||L-^^ , 

||V5,(tH||L- < Ct-i/°||i;o||L- +Ci'"'/"lbo||i, 

M \/ 1 1 rv-\-(i. " " rv-\-n. " n ■ 



With these estimates, we are in a position to construct solutions to 
system (in])-(IIl3]) also in the weighted space L^+d(M'^) C {R'^) n L°° {R'^) . 

Proposition 3.4. Let a € (1,2] and b € (l,d]- Assume that u is the 
solution of system (ll.ip - ()1.3p . constructed in Theorem \2.1\ supplemented 
with the initial condition uq € L'^^^{W^). Then u G C([0, T], L2'+rf(M'^)). In 
particular, for each ^ < a we have J^d {xl"* u{x , t) dx < oo. 

Sketch of proof. Here, it suffices to modify slightly the argument from [1'6\ 
Proof of Prop. 3.3.i] written the case of a convection equation with the frac- 
tional Laplacian. In that reasoning (as well as in the above proof of Theorem 
12. ip . we construct solutions to equation ()3.ip by applying Lemma l3. II with 
the Banach space Xt = C([0, T], L^_^^(M"')). Obviously, y = Sa{-)uo £ Xt 
by Lemma 13.31 Next, we show the following estimate of the bilinear form 
from (13. 6p 

\\H{u,v)\\xr <CT^~^'''\\u\xAv\\xr 

for all u,v £ Xt and a constant C independent of u,v. Here, in order to 
estimate the singular integral operator ()1.3p in the spaces ^^(M'^)), one 
should follow the reasoning from [5l Sec. 2]. Let us omit other details of 
this classical argument. □ 

4 Blow up of solutions 

The main role in our proof of the blow up of solutions to (jl.ip ~ (|1.3p is played 
by the following smooth nonnegative weight function on M'^ 

^{x) = ip^ix) = {1 + - 1 (4.1) 
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with 7 G (1,2]. Since (1 + \x\'^)^ < (1 + l^^j'^)^, we have for each e > 0, 
suitably chosen C{e) > 0, and for every x G M'^ 

ip{x) <\x\^ <s + C{e)ip{x). (4.2) 

Next, let us state two auxiliary results concerning the weight function if 
which will be used in the proof of Theorem l2.3l Here, for a given ip E C^{W^), 
we denote by D'^ip its Hessian matrix. Moreover, the scalar product of 
vectors x, y G M'^ is denoted by x • y. If A is either a vector or a matrix, the 
expression \A\ means its Euclidean norm. 

Lemma 4.1. Let a G (1,2), 7 G (l,a), and if be defined by |-^. j[ ). Then 

(-A)"/V e L°°{^'^). (4.3) 
Proof. First note that by a direct computation we have 

V99(a;) = 7(1 + |xp)2-^x (4.4) 

and 

d,^d,M^) = (7(1 + \x\^)5i, - 7(2 - j)x, xj) (1 + [x|2)i-2. (4.5) 

In particular, for every R > there exists C{R,j) > such that for all 
|x[ > R we have 

|V(^(x)| < C7(i?,7)|xlT-^ and \D^ip{x)\ < C{R,-f)\x\"'~^ . (4.6) 

Now, we apply the following Levy-Khintchine integral representation of 
the fractional Laplacian 

(with a suitable constant C(d, a)) which is valid for every a G (1,2), see, 
e.g. \16\ Th. 1] for a detailed proof of that version of Levy-Khintchine for- 
mula. Using the Taylor expansion and estimates (j4.6p one can immediately 
show that {—A)°'/'^(p{x) is well defined for every x G M'^ and, moreover, 
sup|^l<^ |(— A)°/2(^(2;)| < 00 for each R > 0. 
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In order to obtain an estimate uniform in a; G M"^, we assume that |x| > 1 
and we shall estimate the integral on the right hand side of (I4.7p for \y\ < 
\x\/2 and \y\ > \x\/2, separately. 

If \y\ < I a; 1/2, by the Taylor formula and the second inequality in (14. 6p . 
we obtain 

1 

\ip{x + y)- ip{x) - Vip{x) ■ y\ < -\y\^ \D'^(f{x + sy)\ ds 

< C|yp / \x + syp-^ ds. 
Jo 

Since |y| < \x\/2 and s € [0, 1] we can estimate 

l^^ + syl > ||x| — > ||x| 

Consequently, for 7 — 2 < 0, we obtain 

ip{x + y)- tp{x) - Vip{x) ■ y 



1, , 

> 



l<l^l/2 



< C\x\^-^ 



■dy 



dy 



(4.8) 



\d+a-2 



C\x 



\-f—a 



\y\<\x\/2 \y\' 

for all |x| > 1 and a constant C > independent of x. 

If \y\ > |2;|/2 and |x[ > 1, we combine first inequality from (14. 6p (remem- 
ber that 7 — 1 > 0) with the Taylor expansion to show 

\ipix + y)- ipix)\ < \y\ \V^{x + sy)\ds < C\y\ {Ix^^ + ly^^) . 

Jo 

Hence, 

ip{x + y) - (p{x) - Vip{x) ■ y 



y\>\x\/2 

< cIlxP-^ 



dy 



y\>\x\/2 



+ 



■dy 



l?^l>l^-|/2 



dy 



(4.9) 



\y\d+a- 



7—0 



for all |x| > 1 and a constant C > independent of x. 

Finally, inequalities ()4.8p and ()4.9p complete the proof because 7 < a. 

□ 
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Remark 4.2. Note that above we have, in fact, proved that 

sup (l + l^l""''') (— A)"/^(^(x) < oo for every 7 G (1, a). 



Lemma 4.3. For every 7 E (1,2], the function ip defined in (^TTp is locally 
uniformly convex on Mf^. Moreover, there exists K = K{'j) such that the 
following inequality 

(V^(x) - V^(y)) ■{x-y)> ^ ^ I^IZ~f\y\2-, (4-10) 
holds true for all x, y £ M'^. 

Proof. Using the exphcit expression for the Hessian matrix of ip in (j4.5p we 
obtain 

7(1 + |x|2)|y|2 - 7(2 - 7) (Zixhl + Ei^jXi^jyiVj) 

D"^(-)yy = (1 + 1,|2)2-,/2 (4-11) 

for every x,y £ W^. Now, by the elementary inequahty XiXjyiyj < + 
we immediately obtain 



XiXjViVj < 



2 2 



Consequently, 



Since 7 G (1,2], applying estimate (I4.12|) to (14. lip we get the inequality 

n2 . ^ . 7(l + (7-l)N^)bP 

^^(-^yy^ (1 + ^2-7/2 

which leads directly to the estimate from below 

^V(x)y-y>^^J^i^. (4.13) 
Finally, it follows from the integration of the second derivative of ip that 
(Vip{x) - Vip{y)^ . ix-y) = j D'^Lp{x + s{y - x)) {x - y) ■ {x - y) ds. 
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d 

—w 
dt 



Hence, using inequality (14.130 and the estimate 

(1 + |x + s{y - < dl + + 

valid for all x, y G M*^, s € [0, 1] and a constant C > independent of x, y, s, 
one can easily complete the proof of Lemma 14.31 □ 

Proof of Theorem \2.3[ We consider the function 

w = w{t) = / ip{x) u{x,t) dx, 

where is defined in ()4.1|) and 1 < 7 < a. Note that, in view of inequalities 
()4.2p . the quantity w is essentially equivalent to the moment of order 7 
of the solution u. Moreover, it satisfies the relation 

'-A)°'^^u{x,t)ip{x)dx - / u{x,t)Bu{x,t) ■\/ip{x)dx 
-A)"/V(a;)'u(x,t) dx 

-ifff (v^W-v^M).(.-,) "'"-''"fav'> d.di, 

2 jRd jRd V y \x- yl"' 

(4.14) 

after using the definition of the form B in (II. 3p and the symmetrization of 
the double integral. This computation resembles the usual proof of blow up 
involving the second moments, cf. [21 [TOl [T5l I14j . 

i) For a = 2 = 7 (hence for (/^(x) = |xp) and for f3 = d, the equality 
(|4.14p can be rewritten as follows 

^w{t) = 2dM - Sd,(3M\ 

Evidently, for M > 2d/sd.f3, this implies the equality wiT) = for some 
r > 0, a contradiction with the global existence of nonnegative solutions. 
Thus, we recover the result in [10\ Prop. 4.1] refined in I14j. 

ii) For 1 < /3 < d and fixed M > 0, we are going to use the following 
simple identity 

= / / u{x,t)u{y,t)dxdy 

jR<i JRd 

n X, t mly, t r ; ; dx dy 

^ ^ ^ ^(l + |x!2-7 + |y[2-7)^ |x - y|- 
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1/p' 



with some u > and 5 > 0. We apply now the Holder inequality with the 
powers p > 1 and p' = chosen so that 

vp = d- 13, 5p=l, vp' + (2 - 7) V = 7. (4.15) 

Of course, such a choice of v, 5, p is possible whenever (3 < d and 7 < 2 
because we only need d — /3 + 2 — 7 = j{p — 1). If (3 = d, it suffices to take 
v = and p = 2/7>l. As a consequence, we get 

x([ [ u{x,t)u{y,t)\x-y\''P' {l + \x\^-^ + \yf-^Y^' dxdy) 

(4.16) 

where the integral J{t) satisfies 

^ r r u{x,t)u{y,t) dxdy 

J^J^. \x-y\'i-P l + |x|2-7 + |y|2-7 

<U [ (v,(.)-v^(,)).(.-,)f^^2lM^dxd, 

by Lemma 14.31 

It follows from relations (|4.15|) and inequalities (|4.2|) that there exists 
a constant Ci > such that 

\x - yr' (1 + \x\^-^ + |yp-^)'^' < Ci (1 + (^(x) + ^(y)) . (4.18) 

Hence, ()4.16p implies 

< Cl^^'j{t)^/'P(M^ + 2Mw{t))y^^ . (4.19) 



Going back to identity ()4.14p we obtain from Lemma l4.1l and from inequal- 
ities (fil71) - (jil9]) that 

d 

—wit) < C2M - C3—: — — T (4.20) 

dt + 2Mw{t))p/p' ^ ' 

with C2 = ||(— A)°/^99||oo and a suitable constant C3 > 0. 

Now, we fix for a while M = Mq in (j4.20p so large in order to have 
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Hence, there exists C4 = C4(Mo) > such that for < w{0) < C4 we still 
have 

It is clear that if initially < w{0) < C4 then, by inequality (j4.20p with 
M = Mo, the function w{t) is decreasing in time. Moreover, 

—w(t) < C2M0 - Cs—^ ° , , , < 

dt + 2Mow{0))P/p' 

and, consequently, 'w{T) = for some < T < 00. This contradicts the 
global in time existence of regular nonnegative solutions of ()l.ip - p.3p . Fi- 
nally, note that due to the first inequality in ()4.2I) . it suffices to assume 



w{0) < 



/ |x|'^no(x)dx < C4 and / no(x)dx = Mo, (4.22) 

jRd. jRd 

in order to obtain the blow up in a finite time of the corresponding solution. 

Now, assume that J^iu{x,t)dx = J^aUo{x)dx = M / Mq. Recall 
that system ()l.ip - ()1.3p is invariant under the rescaling (j2.ip . Choosing 
^a+f3-2-d ^ Mq/M we obtain J^aU^ix,t)dx = J^aU^{x)dx = Mq and, by 
(j4.22p . the blow up of the solution takes place under the assumption 

x\'^Uq{x) dx < C4. 

Changing the variables and using the explicit form of A we obtain the blow 
up of solutions to (ll.l|) - ()1.3p under the following assumption on the initial 
condition 

1+: 



-1+ 



d+2-a-/3 



xP'Uo(x)dx < C4M0 «+^-^-^ / uoix)dx 



□ 



Proof of Corollary \2.4\ We follow the proof of Theorem 12. 3[ In particular, 
we choose so large that the inequality (|4.2ip holds true for all Mq > M^. 
This leads to the blow up of the corresponding solution under the assumption 
(j4.22p imposed on the initial data. 
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To complete the proof, we use the scahng argument again. By (12.11 
u^{x,t) = X'^u{Xx, X'^t) is a solution for every A > 0. Note now that 

„5(.)d.= / „„(:.Od.a„d / |.r„J(.)d. = A-/ |.r,.„(.)dx. 



Hence, each initial data uq £ L^{E?, (1 + \x\'~') dx) satisfying Jj^a ^0(2;) dx = 
Mq > leads to the blow up in a finite time of the corresponding solution 
because the moment condition in (|4.22p can be satisfied replacing u by 
and choosing A large enough. □ 
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